function dP = Pint2(y,H,hunt,impossible,phi,s,Vb,Nu)
% Calculate integrand of choice probability 

e = (y + repmat(Vb(1,:,:),[H+1,1,1]) - Vb)./repmat(phi,[1,1,Nu]);
CDF = exp(-exp(-e(1:H,:,:)));
CDF(repmat(impossible(1:H,:),[1,1,Nu])) = 1;

if hunt < H + 1 
    
    CDF(H+1,:,:) = 1 - exp(-s.*e(H+1,:,:));
    dP = prod(CDF).*exp(-e(hunt,:,:)).*1./phi(hunt);
    
else
    
    dP = prod(CDF).*s.*exp(-s.*e(hunt,:,:));
    
end